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We present our preliminary analysis for the chamonium and D s spectra obtained from = 2 dynamical 
anisotropic lattices. We use 12 3 X 80 lattices with lattice spacing at = 7.35 GeV - 1 and anisotropy of six. Meson 
correlators are computed using all- to-all propagators together with variational analysis. 



1. INTRODUCTION 



In recent years, there has been renewed interest 
in charm physics. Many new states such as the 
X(3872) the Y(4260) and the D sJ have been observed 
P, 0, H, 0, and their precise measurement including 
J PC numbers, has become an important topic both 
experimentally and theoretically. 

In principle, lattice QCD should be able to answer 
these questions from first principles but it requires 
high precision numerical simulations. In this region 
of the quark mass relativistic effects could be im- 
portant. However, simulations of the charm quark 
with isotropic lattices are expensive. In this study we 
use a relativistic anisotropic lattice formulation with 
Nf = 2 dynamical quarks to study charmonium and 
D s spectra. In this formulation, the lattice is dis- 
cretized along the spatial, a s , and temporal, at, direc- 
tions with £ = a s /a t 3> 1. Anisotropic lattices have 
the advantage of having small discretization errors in 
the temporal direction whilst keeping the computa- 
tional cost down. In addition, keeping a small tem- 
poral lattice spacing allows us to increase the number 
of time slices which in turn makes it easier to iden- 
tify the plateau regions in effective mass plots. It is 
difficult to achieve this using isotropic lattices since 
the heavy hadron correlators with a charm or bottom 
quark have a signal that decays rapidly. Our aim is 
to be able to extract the excited spectra with small 
statistical errors. 

Our lattice actions are described in section [5] In 
section [3] we summarize the parameters and the inter- 
polating operators used in the simulation. In section^ 
we discuss the method to obtain excited spectra. In 
conclusion, we discuss our results and future plans. 



2. ANISOTROPIC ACTIONS 

In this section, we describe the gauge and quark 
actions used in the simulation. The gauge action is 
a Two-Plaquette Symanzik-improved action which is 
designed to study glueballs @|. It is given by 
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where f2 s and O t are the spatial and temporal 



plaquettes respectively. fl 



and fij^ refer to 



space-space and space-time rectangles and fii 2 ^ — 
I J2x i>j [l ~ Pij { x )Pij ( x + t)\ ■ This action has lead- 
ing discretization errors of 0(cQ, a 2 , a s a 2 s ). 

The 0(al, a t , a s a s )-improved anisotropic quark ac- 
tion 0] is given by 
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where r and s are the Wilson parameters and fi r = 
(1 + ratTOo/2). The derivatives are defined as 

V^VO) = [U^tpix+fi) - Ufa-jiMx-fi)] , 
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and £° is the bare quark anisotropy. In order to max- 
imize the plaquette stout links Q are used. We used 
the same quark action to simulate light sea quarks 
and heavy valence quarks. Our sea quark mass in this 
simulation is close to the strange quark mass. 
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Table I Simulation parameters. 



Configurations 


250 (a t m c = 0.117, 




a t rn sea = atmu g ht = —0.057) 


Dilution 


time, space(even/odd), color (cc) 




time, space(even/odd) , (D s ) 


Physics 


S, P and D waves, hybrids 


Volume 


12 3 x 80 


N f 


2 


a s 


~ 0.17 fm 


ar 1 


7.35 ± 0.03 GeV 



The ratio of the lattice spacings, a s /a t , which ap- 
pears in both the gauge and the quark actions, as 
£® and £° respectively, are bare parameters that need 
to be tuned. In a quenched study, this tuning can be 
done separately. However, sea quark effects in dynam- 
ical simulations lead to a simultaneous non-trivial tun- 
ing of the anisotropies. This procedure is explained in 
detail in Ref. 9]. For the results presented in these 
proceeding the renormalized anisotropy, £r is set to be 
6. 



3. SIMULATION DETAILS 

In this study, we obtained the charmonium and the 
D s spectra from 12 3 x 80 lattices with 250 configu- 
rations. We tuned the valence charm quark mass to 
atm c = 0.117 in order to get the J/ip mass correct 
while the light quark mass is a t m sea = at m light = 
—0.057 which is close to the strange quark mass. We 
use the all-to-all propagator method with "dilution" of 
Ref. fioj with no eigenvalues for the charm quark prop- 
agators and 20 eigenvalues for the strange quark prop- 
agators. We used time, space(even/odd) and color di- 
lution for the charmonium study while color dilution 
is omitted in the D s case. The lattice spacing is set 
from the spin-averaged (IP-IS) splitting in the char- 
monium system and found to be at — 0.0272 fm. The 
parameters used are listed in table |TJ 

Taking advantage of the all-to-all propagators, we 
use a variational method [ll], E3] m order to get a bet- 
ter overlap with higher excited states where we use a 
spatially extended operator basis [TH . The lattice op- 
erators used in this study along with their putative 
continuum spin assignments are given in table [TTJ We 
have assumed the simplest possible continuum assign- 
ment is correct. The validity of this assumption is 
under investigation. Using the operator basis in ta- 
ble [HI for a given channel R one can construct the 
matrix 

Cif(t) = (0\Oi R \t)Ol^(0)\0), (5) 



where a, (3 — l,...,n represent the different inter- 
polating operators, constructed by applying different 
levels of quark smearing. The different energy levels 
can then be obtained from 

lim X a (t, t ) = e -^- u ^ \l + 0(e- tAE *)} , (6) 

where X a are the eigenvalues of the matrix 
C(^o) -1 ^ 2 C(£)C(£o) -1 / 2 and i is some small reference 
time. We performed single state fits to the diagonal 
elements in order to extract the ground and excited 
states. 



Table II The operator basis used to obtain the cc and D s 
spectra. Definitions of the Si,pt and t% are given in 
Ref. 0. 





2S+l Lj 


STATE 


OPERATORS 


0"+ 

\— 


'So 
3 Si 


Vc,Vc 

JM^(2S) 


75,75 J2i s » 

7i,7j J2i s i 


1+- 

0++ 

1++ 

2++ 


'Pi 

3 Po 

3 P1 
3 P 2 


h c ,h c 

XCQ tXcQ 


1,7 -p 

7s7i,7 x P 
7 x p , 71P1— 72P2 
273P3 — 71P1 — 72P2 


2- + 
2"" 

3— 


'D 2 
3 D 2 

3 D 3 


l'D 2 
l 3 D 2 

1 3 D 3 


7s(«i - s 2 ),75(2s 3 — si — s 2 ) 
7j(si - s fc ),7iti - 72*2, 
273*3 - 71*1 - 72*2 
7-f 


1-+ 


Hybrid 


qqg 


7 x u 



4. ANALYSIS 

Time diluted all-to-all propagators introduce ran- 
dom noise at each time slice which makes it difficult to 
indentify a plateau region in the effective mass plots 
which fluctuate more than point-to-all propagators. 
However, we fit the correlators and a better picture 
can be obtained from "sliding window" plots (or t m i n 
plots). For a fixed value of some t max , we vary the 
tmin value and plot the fitted mass. This is illustrated 
in Figs. [T]and [2]for the J ftp and D s - and 1~ states, 
respectively [Tilflll . We choose our fits based on the 
X 2 /Ndf (<2), fit range (where the fits are stable) and 
the fit quality (Q>0.2). These values are chosen from 
our earlier simluations with smaller lattices. Most of 
our results have better \ 2 and Q values. 

We performed the same analysis for the cc and D s 
systems. Our preliminary spectra are shown in fig- 
ures [3] and [U 
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Figure 1: Sliding window plots for the J/ip obtained 
from the variational analysis. The plot shows the ground 
state in the lattice Tf irreducible representation which 
corresponds to the continuum J = 1 state. The higher 
lying states determined in this channel are assumed to be 
radial excitations of the J/$. This is under investigation. 
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Figure 2: Sliding window plot for the D s and 1 
states. The corresponding experimental value of the 
1" - 0" splitting is 143.8MeV. 



5. Conclusions and Outlook 

We have presented our preliminary results for the 
charmonium and D s systems from Nf — 2 dynamical 
anisotropic lattices. All-to-all propagators are essen- 
tial in this study and allow us to use a wide range of 
operators and the variational analysis. For the char- 
monium system we have good signals for the S, P and 
D waves and the 1 ^ hybrid. We are planning to ex- 
pand the simulation to include the 1 D- waves and 
other hybrids. We found the hyperfine splitting in this 
system to be small. The effect of the chromomagnetic 
term, cbX • B, disconnected diagrams and stout link 
smearing are being investigated as possible reasons for 
this. The D s system is simulated with a low level of di- 
lution. A simulation with a higher level of dilution for 
this system, with a wider range of operators, is being 



Figure 3: Preliminary charmonium spectrum for the S, 
P, D waves and the hybrid 1 . The results of this study 
are the blue and red bands. The highest lying radial 
excitations identified in each channel are coloured red to 
indicate that these are not free of further excited-state 
contamination. 
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Figure 4: Preliminary D s spectrum for the S and P 
waves. The results of this study are shown by the blue 
and red bands. As in the charmonium case the highest 
lying radial excitations identified in each channel are 
coloured red to indicate that these are not free of further 
excited-state contamination. The blue dots represent the 
UKQCD (N f = 2) results 0. 



performed. Both simulations are performed at single 
lattice spacing where the sea quark mass is around 
the strange quark mass. Simulations with finer lat- 
tices spacings are currently under investigation. The 
results clearly demonstrate the power of the all-to-all 
propagators combined with a variational analysis. We 
have extracted a large number of orbital and radial 
excitations in the charmonium and D s systems. Fur- 
ther work is also underway to address the continuum 
spin-identification of these lattice determinations. 
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